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^ ■ Abstract 

f~^ ' We discuss new diffr active mechanism of central exclusive production of W~^W pairs in proton- 

CN . proton collisions at the LHC. We include diagrams with intermediate virtual Higgs boson as well as 

F^ . quark box diagrams. Several observables related to this process are calculated. Predictions for the 

[t \ total cross section and differential distributions in T^-boson rapidity and transverse momentum as 

well as WW invariant mass are presented. We also show results for different polarization states of 
the final W bosons. We compare the contribution of the 77 — )■ W~^W~ mechanism considered 
in the literature with the contribution of the diffractive mechanism through the gg — >• W^~''VF~ 
i-C . subprocess for the different observables. The phase space integrated diffractive contribution when 

I \ separated is only a small fraction of fb compared to 115.4 fb of the 77-contribution without ab- 

rr^l sorption. The latter contribution dominates at small four-momentum transfers squared in the 

,^ ' proton lines and in a broad range of W^W~ invariant masses. This offers a possibility of efficient 

searches for anomalous triple-boson (714/^VF) and quartic-boson (77VFVF) couplings and testing 
CN , models beyond the Standard Model. We discuss shortly also the pp — ;■ OT77 process, where the 

p^ I box contribution is very similar to that for W^W~ and compare our results with recent CDF data. 

cn ' Nice agreement has been achieved without additional free parameters. 
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I. INTRODUCTION 

The central exclusive production (CEP) process pp ^ p + X + p, where X stands for a 
centrally produced system separated from the two very forward protons by large rapidity 
gaps, has been proposed in Refs. [l|, \^ as an alternative way of searching for the neutral 
Higgs boson (see Ref. |3| for a review). If momenta of the outgoing protons are measured 
by forward proton detectors placed at 220 m and 420 m from the ATLAS/CMS interaction 
point 1^], the mass of the X system may be reconstructed [Si] with very precise resolution. 

The exclusive reaction pp — )■ pHp has been intensively studied by the Durham group 

in the last decade. This study was motivated by the clean environment and largely 
reduced background due to a suppression of bb production as a consequence of the spin-parity 
conservation in the forward limit. However, very recent precise calculations of Refs. [Tj] have 
shown that the situation with Higgs CEP background in the hh channel is more complicated 
and the signal is to a large extent shadowed by the exclusive non-reducible continuum bb 
production. In addition, reducible backgrounds from a misidentification of gluonic jets as 
6-quark jets can be very difficult to separate [§]. Since the total cross section for the Higgs 
CEP is quite small and rather uncertain, the issue with the Higgs CEP is still far from its 
final resolution, from both theoretical and experimental point of view. 

The final system X in the midrapidity region is predominantly produced in the Jz = 
state as dictated by the well-known J^ = selection rule |6|. However, corrections to this 
rule due to slightly off-forward protons can be important for lower (a few GeV) mass central 
systems and may lead to sizeable contributions in the observable signals, in particular, in 
the Xc mesons [a, llO|, hh [Tj and gg [g, [ll| CEP. The emission of gluons from the "screening" 



gluon could also violate the J^ = selection rule as has recently been emphasized in Ref. |ll| . 

In order to reduce the theoretical uncertainties of the CEP mechanism, coming from both 

the hard subprocess (Sudakov form factor Uj, [l^, next-to-leading order QCD corrections 



[13(1) and the soft interactions (color screening effects at extremely small gluon x [1^, rapidity 
gap survival factor [l5|, poorly known unintegrated gluon distribution functions (UGDFs) 
at small gluon g^ and x |9|, [lO|), new experimental data on various exclusive production 



channels are certainly required and expected to come soon from ongoing LHC measurements. 
In particular, as it was stressed e.g. in Ref. [l2| the measurements of the exclusive dijets 
production at the LHC could largely reduce the theoretical uncertainty in the Higgs boson 
CEP. Other measurements, e.g. heavy quarkonia |9|, llO|, 77 10|, high-p^ light mesons 



16|, |17[, associated charged Higgs H^W [18|] CEP, etc., are also important in this context. 



Some of these results have been compared to experimental data from the Tevatron [19|, |20 
and a rough quantitative agreement between them has been achieved. 

In this paper, we focus on exclusive production of W^W~ pairs in high-energy proton- 



proton collisions. It was found recently [21|, |22| that the reaction is an ideal case to study 



experimentally '~fW~^W~ and '-f'-fW~^W~ couplings ^. 

The 77 — )■ W~^W~ process is interesting reaction to test the Standard Model and any 
other theory beyond the Standard Model. The linear collider would be a good option to 



study the couplings of gauge bosons in the distant future. For instance in Ref.[2J] the 
anomalous coupling in locally SU(2) x U(l) invariant effective Lagrangian was studied. 
Other models also lead to anomalous gauge boson coupling. 

The photon-photon contribution for the purely exclusive production of W^W' was con- 



^ Some more subtle aspects of the beyond Standard Model anomalous couplings were discussed e.g. in 23|. 



sidered so far in the literature. The diffractive production and decay of Higgs boson into 



the W^W~ pair was discussed in Ref. J25|, and the corresponding cross section turned out 
to be significantly smaller than that for the 77-contribution. Provided this is the case, the 
W~^W~ pair production signal would be particularly sensitive to New Physics contributions 



in the 77 — )■ W^W~ subprocess |2ll . |22| . Similar analysis has been considered recently for 
77 — )■ ZZ |26|. These previous analyses strongly motivate our present detailed study on a 
competitive diffractive contribution. The pp — )• pW^W~p process going through the diffrac- 
tive QCD mechanism with the gg — )■ W^W^ subprocess naturally constitutes a background 
for the exclusive electromagnetic pp — ?■ ^(77 — t- W'^W~)p process. We consider not only 
the mechanism with intermediate Higgs boson but also quark box contributions never esti- 
mated in exclusive processes. Both the Higgs and box contribution may interfere together. 
We discuss here the interference effects. Corresponding measurements will be possible to 
perform at the ATLAS detector with the use of very forward proton detectors 2l| . In order 
to quantify to what extent the QCD mechanism competes with the "signal" from the 77 
fusion, we calculate both contributions and compare them differentially as a function of 
several relevant kinematical variables. 

Since the box contribution of exclusive diffractive pp — )• ppW^W~ process is very similar 
to the pp — )■ PP77 process which has been measured recently |20| , we discuss the latter one 
and compare corresponding results with the recent CDF data. 

II. DIFFRACTIVE MECHANISM OF EXCLUSIVE W+W PAIR PRODUCTION 



A schematic diagram for central exclusive production of W^W^ pairs in proton-proton 
scattering pp — )■ pW^W^p is shown in Fig. [1] Similar mechanisms have been considered in 
inclusive production of W^W~ pairs (see e.g. Refs. |27|, |29|, |30| ) . In what follows, we use the 
standard theoretical description of CEP processes developed by Khoze, Martin and Ryskin 
for the exclusive production of Higgs boson 




FIG. 1: Generic diagram for the central exclusive WW pair production in pp collisions. Momenta 
of incident particles are shown explicitly. 

The momenta of intermediate gluons are given by Sudakov decompositions in terms of 
the incoming proton four-momenta pi,2 



qi = xipi + gi_L, q2 = X2P2 + q2±, < a;i,2 < 1, 
qQ = x'pi-x'p2 + qo±-qo±, a;' < Xi,2, 



(2.1) 



where xi,2, x' are the longitudinal momentum fractions for active (fusing) and color screening 
gluons, respectively. 



In the forward proton scattering limit, we have 

tl,2 = {Pl,2 - Pl,2? - P'l,2± -^ , 

q± = qo± - -qi± = q2± ■ (2.2) 

The QCD factorisation of the process at the hard scale fxp is provided by the large invariant 
mass of the WW pair Mww, i-e. 

fip = s X1X2 ~ M^yy . (2.3) 

It is convenient to introduce the Sudakov expansion for W"^ boson momenta 

k+ = xfpi+ X2P2 + k+±, k_ = x^pi + X2P2 + k__L (2.4) 

leading to 

^1,2 = ^12 + ^1,2' -^1,2 ~ 7^^ ' "^1.2 ~ 7^^ ' ^^±± ~ '^^W "I" l^'i^l 1 (2.5) 

in terms of W^ rapidities y± and transverse masses m±_L. For simplicity, in actual calcula- 
tions we work in the forward limit given by Eq. (12. 2p . which implies that fc_|__|_ = — fc__|_. 

In actual calculations below, W^ bosons are assumed to be on-mass-shell, whereas par- 
ticular contributions to the observables can then be estimated in the narrow-width approxi- 
mation. For example, in the leptonic channel we have the following observable cross section 

ai+^i-^ - aww X BR{W+ -^ l+v) BR{W- -^ I'v) , (2.6) 



where BR(W^+ -)■ l^v) = (10.80 ±0.09) x 10~^ [lH for a given lepton flavor. Both electrons 
and muons can be used in practice |2ll |. 



We write the amplitude of the diffractive process, which at high energy is dominated by 
its imaginary part, as 

KA / J. J. \ • ^ /" j2 T/ / ; ; N /g(?0;?i;^l)/g(Q'0;?2; ^2) /^ „x 

Mx+x^{s,ti,t2) ~ ts— / d'qo±Vx^x-{qi,q2,k+,k_)^ 2 tS ' (2-7) 

^ J 9o± Ql± Q2± 

where X± = ±1, are the polarisation states of the produced W"^ bosons, respectively, 
fg{.'^i^'^2]'t) is the off-diagonal unintegrated gluon distribution function (UGDF), which de- 
pends on the longitudinal and transverse components of both gluons momenta. The gauge- 
invariant gg — > W^ W^_ hard subprocess amplitude Va+a_('?i, ^2, fc+, fc-) is given by the 
light cone projection 

Va, A_ = n n Vx,X^,au = Vx,X.,av, qiVx+X.,av = Q'2^A+A_,ui/ = 0, (2.8) 

'^ S Xi X2 

where n^ = Pi^2/Ep,cms and the center-of-mass proton energy Ep^cms = V^f^- We adopt 
the definition of gluon transverse polarisation vectors proportional to the transverse gluon 
momenta qi,2±, i-e. ei,2 ~ qi,2±/xi^2- The helicity matrix element in the previous expression 
reads 

V,^X{quq2,k^,k_) = e*'''(A;+, A+)e*''^(fc_, A_)IX% (2-9) 



in terms of the Lorentz and gauge invariant 2 — )■ 2 amplitude V^J:' and W boson polarisation 
vectors e{k, A). Below we will analyze the exclusive production with polarized W~^ or W~. In 
Eq. fl2.9p e^(/i;+, A+) and ei,(fc_, A_) can be defined easily in the proton-proton center-of-mass 
frame as 

e(fc,0) = I — — , cos(f)sm9, sm(f)sm9, cosO 

e{k,±l) = —= {0, isincf)^ cos9cos(f), —icos(f)^ cos9sm(f), ±sm9) , (2.10) 

V2 

where is the azimuthal angle of a produced boson, and satisfy e^(A)e*(A) = —1 and 
e*(fc+, X+)k1 = el{k_, XJ)k^_ = 0. In the forward limit, provided by Eq. (12. 2p . the azimuthal 
angles of the W^ and W~ bosons are related as 0_ = 0+ + vr. 

The diffractive amplitude given by Eq. (12. 7p is averaged over the color indices and over 
the two transverse polarizations of the incoming gluons. The relevant color factor which 
includes summing over colors of quarks in the loop (triangle or box) and averaging over 
fusing gluon colors (according to the definition of unintegrated gluon distribution function) 
is the same as in the previously studied Higgs CEP (for more details on derivation of the 
generic pp — )■ pXp amplitude, see e.g. Ref. [Sj). The matrix element Va+,a_ contains twice 
the strong coupling constant g"^ = Anag- In our calculation here we take the running 
coupling constant asifJ^lard — ^ww) which depends on the invariant mass of WW pair as a 
hard renormalisation scale of the process. The choice of the scale approximately introduces 
roughly a factor of two model uncertainties when varying the hard scale ^hard between 
2Mww and My^w /2 values. 

The bare amplitude above is subjected to absorption corrections that depend on the 
collision energy and typical proton transverse momenta. As in the original KMR calculations 
[61], the bare production cross section is usually multiplied by a rapidity gap survival factor 
which we take the same as for the Higgs boson and hh production to be Sg = 0.03 at the 
LHC energy (see e.g. Ref. 



A. The hard subprocess 

The typical contributions to the gg — )■ W^W~ subprocess are shown in Fig. [2l The total 
number of topologically different loop diagrams amounts to two triangles, and six boxes. 
In the central exclusive W^W~ production, triangle diagrams with 7 and Z bosons in the 
intermediate state are suppressed due to the J^ = and parity selection rule for singlet 
gluon-gluon to (virtual) photon transition strictly valid in the on-shell limit of fusing gluons 
and Landau- Yang theorem for intermediate Z boson. 

Then the only non-zeroth contribution comes from the Higgs resonant diagram, and in 
the next subsection we will discuss it in detail. However, this can only lead to a sizeable 
enhancement of the cross section close to its threshold rriho — Mww ^ 2mw 25|. The 



Standard Model Higgs bosons with such large masses have been recently excluded by the 



Tevatron [32| and LHC [33|, |3J] measurements. For yet allowed values of Higgs mass 115 
GeV < rrtho < 130 GeV, corresponding contribution to the W^W~ channel is far from the 
Higgs boson resonance and turned out to be suppressed compared to box contributions at 
low invariant masses. However, due to interference effects at rather large invariant masses 
M\Yy/ the resonant (triangles) contribution could become comparable to the non-resonant 



(boxes) one. Below, for comparison we have calculated box and triangle (through the s- 
channel SM Higgs boson exchange) contributions in different phase space regions ^ which 
could be interesting for future measurements with forward detectors at ATLAS or CMS. 
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FIG. 2: Representative diagrams of the hard subprocess gg — ?■ VF VF^, which contribute to the 
exclusive WW pair production. 



1. Higgs contribution 

The matrix element for the gg ^ h^ ^ W~^W~ transition with intermediate s-channel 
Higgs boson exchange (see first two diagrams in Fig. [2]) can be written in the narrow-width 
approximation as 



Vgg^h0-,w+w-iqi,q2,k+,k-) = 5^^\qi + q2 - k+ - k. 



X 



ygg^h'i{(il-,<l2-,Ph('] 



Mir,,r - mi, + iM^ 



WWi- tot 



:;j;-VhO^W+W-{k+,k_,X+,X_), (2.11) 



where the Higgs boson momentum is Pho = qi + q2, and the 5- function reflects the momentum 
conservation in the process. In order to get a correct resonant invariant mass distribution, 
the standard Breit-Wigner Higgs propagator with the total Higgs decay width F(*q^, which 



can be found e.g. in Ref. [36j, is used. 

In Eq. (12. lip , first the gg — >■ h^ amplitude of the Higgs boson production through the 
top-quark triangle in the /ci-factorisation approach can be written as (see e.g. Ref. [37|) 



K 



gg^hO 






"top 



Gf^ 



-1/2 



The second tree- level h^ -^ W^W "decay" amplitude reads: 

e 



V^ 



hO^W+W- 



iniw 



sin^^y 



e*(A;+,A+)e*(A;_,A_) 



(2.12) 



(2.13) 



Close to the VFVF-threshold instability of W bosons [35.] should be included. 



where the polarisation vectors in the direction of motion of W~^ and W~ bosons in the 
proton-proton center-of-mass frame are used in practical calculations. 

Potentially interesting contribution could come from the Higgs resonance if the Higgs mass 
was close to the WW production threshold. Similar resonance effects have been considered 
recently in inclusive |38| and exclusive associated [.18| charged Higgs boson production, and 
large contributions beyond the Standard Model were found. However, the SM Higgs mass 
~160 GeV has been recently excluded in inclusive searches by the CDF Collaboration at 
Tevatron J32l] and by the ATLAS and CMS Collaborations at LHC |33l.l34|. so yet realistic SM 



Higgs boson mass interval rriho ~ 115 — 130 GeV leads to a suppressed triangles' contribution 
to exclusive W^W~ pair production. In the calculation presented here we take niho = 120 
GeV. Since the Higgs mass is certainly much smaller than the threshold value a precise 
value of the Higgs boson mass is not very important. A contribution from an extended 



Higgs sector beyond the Standard Model [38j could be interesting, but we postpone this 
issue for a later study. 

In this work, we are primarily interested in estimation of dominant box contributions as 
well as in possible box-triangle interference effects within the Standard Model as potentially 
important irreducible background for the 77 — >■ W^W~ signal relevant for a precision study 
of anomalous couplings. Thus, our numerical estimates provide minimal limit for the central 
exclusive WW production signal. 

2. Contribution of box diagrams 

The box contributions to the gg — )■ W^W~ parton level subprocess amplitude (see dia- 
grams No. (3-8) in Fig. [2]) for on-shell fusing gluons were calculated analytically by using the 
Mathematica-based FormCalc (FC) |39j package. The complete matrix element was gener- 
ated automatically by the FC tools in terms of one-loop Passarino-Veltman two-, three- and 
four-point functions and other internally-defined functions (e.g. gluon and vector bosons 
polarisation vectors) and kinematical variables. 

At the next step, the Fortran code for the matrix element was generated, and then used as 
an external subroutine in our numerical calculations together with other FC routines setting 
up the Standard Model parameters, coupling constants and kinematics. Instead of built-in 
FC polarisation vectors we have used transverse gluon polarisation vectors which enter the 
projection in Eq. (12. Sp . and the standard W^ polarisation vectors defined in Eq. (12.101) . 
giving us an access to individual polarisation states of the W bosons. In accordance with 
the fcf-factorisation technique, the gauge invariance of the resulting amplitudes for the on- 
mass-shell initial gluons is ensured by a projection onto the gluon transverse polarisation 
vectors proportional to the transverse gluon momenta gi,2± according to Eq. (12. 8p . 

For the evaluation of the scalar master tree- and four-point integrals in the gluon-gluon 
fusion subprocess we have used the LoopTools library [39]. The result is summed up over all 
possible quark flavors in loops and over distinct loop topologies. We have also checked that 
the sum of relevant diagrams is explicitly finite and obeys correct asymptotical properties 
and energy dependence. It is worth to mention that a large cancelation between separate 
box contributions in the total sum of diagrams takes place, which is expected from the 
general Standard Model symmetry principles ^. 



We are thankful to Prof. O. Nachtmann for an enlightening discussion on this matter. 



As soon as the hard subprocess matrix element (denoted above as Va+a_) has been defined 
as a function of relevant kinematical variables (four-momenta of incoming/outgoing parti- 
cles), the loop integration over go± in Eq- (12. 7p was performed to obtain the diffractive am- 
plitude, which then has been used to calculate the differential distributions for (un)polarised 
W bosons. 

As we will demonstrate below, in the Standard Model the total box contribution is 
somewhat larger than the triangle one, for the realistic Higgs boson masses. We, however, 
keep both the triangle and box contributions and investigate a possible interference between 
them, which, in fact, is quite important, especially at rather large iy"'"iy~-pair invariant 
masses, i.e. in the region we are interested in. 

B. Exclusive pp — )• pp77 process 

The same formalism as described above is used to calculate the amplitude for the pp — )■ 
PPIl process. We write the amplitude of the diffractive pp — )■ pp'j'j process as 

Mx+x.{sMM) ^ ^s— / rf Qox^A+A- (?i>?2,A;+,A;_)^^ 2 T^ .(2-14) 

^ J %± Ql± Q2± 

where now X± = ±1 are the helicity polarisation states of the produced photons and corre- 
sponding polarisation vectors are defined easily in the pp center-of-mass frame 

e(A;, ±1) = ^p (0, isin0 =F cos6'cos0, — icos0 =F cos6'sin0, ±sin6') . (2.15) 

v2 

The typical contributions to the leading order gg — )■ 77 subprocess are shown in Fig. |3l 
The total number of topologically different loop diagrams in the Standard Model amounts 
to twelve boxes. So the 77 does not exhibit resonant features, and can potentially serve as 
a probe for New Physics resonant contributions. 
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FIG. 3: Representative diagrams of the hard subprocess gg — >• 77, which contribute to the exclusive 
77 pair production. 

The box contributions to the gg -^ 77 parton level subprocess amphtude in Fig. |3]for on- 
shell fusing gluons were calculated analytically by using the Mathematica-based FormCalc 
(FC) |39| package. The complete matrix element was automatically generated by FC tools 
in terms of one-loop Passarino-Veltman two-, three- and four-point functions and other 
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gluon and vector bosons polarisation vectors) and kine- 

produc- 



internally-defined functions (e.g 
matical variables. 

Other details of the calculation are very much the same as those for the W^W 
tion. We will not repeat here the details. 



C. Gluon fc^-dependent densities in the forward limit 



In the fci-factorisation approach, the density of gluons in the proton is described in terms 
of the off-diagonal unintegrated gluon distribution functions (UGDFs) /g((?o; Q'1,2; ^1,2) = 
fg^{x', a;i_2, Qo_L, 9i/2±; /^fj ^1,2) at the factorization scale fip ~ M^y^r ^ \qo±\- In the forward 
scattering (see Eq. (12.21) ) and asymmetric limit of x' ^ xi_2, the off-diagonal UGDF is 
written as a skewedness factor Rg{x') multiplied by the diagonal UGDF, which describes the 
coupling of gluons with longitudinal momentum fractions Xi^2 to the proton (see Refs. 401.141 
for details). The skewedness parameter Rg is expected to be roughly constant at LHC 
energies and gives only a small contribution to the overall normalization uncertainty. We 
take Rg = 1.3 in practical calculations. In the kinematics considered here, the unintegrated 
gluon density can be written in terms of the conventional integrated gluon distribution 



g{x,ql) as 



a 



fg{qo,qi,2]h,2) ^ Rgfg{xi,2,ql,f^l)exp{bti,2/2) 

d 



Rr 



dlnq"] 



x,,29ix,,2,ql)jTg{ql,fil) exp(6ti,2/2) , (2.16) 



Tg is the Sudakov form factor which 



where the diffractive slope is taken to be 6 = 4 GeV^^ 

suppresses real emissions from active gluons during the evolution, so that the rapidity gaps 

are not populated by gluons. It is given by J4l| 



^9(91,/^!) = exp 



^'F dk^, 



«, 



(fcl) 



k' 



2n 



^ zPgg{z) + J2P,ai^) d^, (2.17) 



where A in the upper limit is taken to be 42 



A = 



\k 



'± 



\k±\ + Mww 



(2.18) 



In our calculations we take /i|^ = M^^,. The choice of the scale introduces uncertainties 
roughly of about factor two. Since in the present calculations we need values of Tg{q\, ^y) 
for extremely large scales ix\ the integration in Eq. (I2.17P is performed rather in \ogiQ{k'^ / k"^) , 
where /co = 1 GeV was introduced for convenience. 



III. FOUR-BODY PHASE SPACE IN THE FORWARD LIMIT 



The diffractive WW CEP amplitude (12. 7p described above is used now to calculate the 
corresponding cross section including certain limitations of the phase space. The cross 
section for the two-boson production can be obtained by integration over the four-body 



phase space given by 

(2vr)' fTTTJ^.A, , . . , , ^ ^M d% d^k^ d'k. 



^ =^^J\M\^\pr+p2-p[-p', -k+- k_) 



{27Ty^2E[ (27r)32E^ (27r)32E+ (27r)32E_' 

(3.1) 

± 



where E[ 2 and E± are the energies of the final-state protons and produced W bosons, 

respectively, |A4P = ^^^ a_ -^a+a_A^a+a_ assuming, as usual, that the helicities of both 
protons are unchanged in the considered process. In order to calculate the total cross section 



one has to take the eight- dimensional integral numerically (for details see e.g. Ref. |43|). 
However, the evaluation of the corresponding hard subprocess amplitude Va+a_ , its subse- 
quent convolution with the gluon UPDFs in the diffractive amplitude (12.71) and the full phase 
space integration (13.11) is extremely time consuming. Clearly the calculation of diffractive 
mechanism must be simplified to be feasible. Such a simplification seems possible for the 
diffractive process considered here. We start from the choice of integration variables as in 
Ref. |43i]. Then 

da = — |A^|2 -^-^^-^ dtidt2d(j)id(j)2 -^J~^ dy+dy_dp^±d(p,n , (3.2) 

where pm± = |fc+i. — fc-±| is the difference between transverse momenta of W~^ and W~, 
k^± and fc__L, respectively, and (pm is the corresponding azimuthal angle. For the sake 
of simplicity, assuming an exponential slope of ti/t2-dependence of the KMR UGDFs (see 
Eq. (I2.16P ). and as a consequence of the approximately exponential dependence of the cross 
section on ti and t2 (proportional to exp(6ti) and exp(6t2)), the four-body phase space can 
be calculated as follows 



1 



2s' 



11 1 1 1 /o ^2Pm± 



da^ -\M\^ (2ny ^J-' dy^dy^dp^^dc^^ . (3.3) 



=0 24 (27r)8 E[E'2 4 62 ' ' 4 



Since in this approximation we have assumed no correlations between__outgoing protons 
(which is expected here and is practically true for the production of bb [7] or gg |8| dijets) 
there is no dependence of the integrand in Eq. (13. 3p on 0^? which means that the phase space 
integration can be further reduced to three-dimensional one. The Jacobian J' in Eq. (13. 2p 
is given in Ref. J43| 



J 



Viz P2z 



1^1 + p'lz \hK + P'lz 



(3.4) 



In actual calculations below we shall use the reduced form of the four-body phase space 
Eq. (13. 3p . and it is checked to give correct numerical results against the full phase space 
calculation for some simple reactions. Different representations of the phase space depending 
on a particular kinematical distributions needed can be found in Ref. |43|. 

IV. 77 -^ W+W- MECHANISM 

In this section, we briefly discuss the 77 — )■ W^W~ mechanism, considered already in 
the literature (see Refs. |2l|, l22|). The relevant subprocess diagrams are shown in Fig. |H 
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Let us start from the reminder about the 77 — )■ W^W~ couphng within the Standard 
ModeL The three-boson WW'j and four-boson WW'yy couphngs, which contribute to the 
77 — )• W^W~ process in the leading order read 



C-ww,, = -e\W-W^^A,A^ - W-A^W^A^) , 



(4.1) 
(4.2) 



where the asymmetric derivative has the form X d'^ Y = Xd^Y — Yd^X. 




7 7 

(11 (21 (31 

FIG. 4: The Born diagrams for the 77 — >■ W W^ subprocess. 

Then within the Standard Model, the elementary tree-level cross section for the 77 — >■ 
W^W~ subprocess can be written in the compact form in terms of the Mandelstam variables 
(see e.g. Ref. |30|) ^ 



da 3a2/3 {^ 2s(2s + 3m^) 



dVL 



2s 



2g^(g^ + 3m^) 
3(m^-t)(m^-M) ■ 3{ml, - i)^{ml, - u) 



+ 



(4.3) 



where /3 = a/1 — 4m^/s is the velocity of the W bosons in their center-of-mass frame and 
the electromagnetic fine-structure constant a = e^/(47r) ~ 1/137 for the on-shell photon. 
The total elementary cross section can be obtained by integration of the differential cross 
section above. 

In the Weizsacker- Williams approximation, the total cross section for the pp — )■ pp(77) — >■ 
W^W~ can be written as in the parton model 



a 



/ dxidx2 f^^ (xi) f. 



WW I \ - I '-\ 

2 \^2)(y-t-t^W+W~\S) 



(4.4) 



We take the Weizsacker- Williams equivalent photon fluxes of protons from Ref. |44l |. 

To calculate differential distributions the following parton formula can be conveniently 

used 

da 1 



(WW 



dy+dy^d^PwL 



IGvr^s^ 



Xlfi {,Xi)X2f2 {x2)\M^^^w+W-{.S,t,U 



(4.5) 



This formula does not include the process with virtual Higgs boson 77 — > i/ -^ W^W 45|. For heavy 



Higgs boson, this would lead to clear Higgs boson signal modifying the cross section (typical resonance + 






background effect) [30j, however, with the present limits for Higgs boson mass 33|,|3j] only deeply off-shell 
Higgs boson contribution could be possible. Also, the diagram with an intermediate Higgs boson is, of 
course, of a higher order compared to the contributions considered here. This automatically means rather 
small effect on the measured cross section, in particular, on the W^W~ invariant mass distribution in 
our case of the four-body pp -^ pW^W~p reaction. 
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where momentum fractions of the fusing gluons x\^2 are defined in Eq. (12. 5p . We shall not 
discuss here any approach beyond the Standard Model. A potentially interestin g H iggsless 



scenario of the lyVT-pair production has previously been discussed e.g. in Refs. [2l|, |22 . 

In Fig. Owe show distribution in ^i = \o'g^^{x\) and ^2 = logio(^2) at a/s = 14 TeV. 
We observe a maximum of the cross section at ^1, .^2 ~ —2 which means that corresponding 
longitudinal momentum fractions carried by photons are typically 10~^. 




00 

^ -1- 




- pp ^ pp W^ W" (via yy) 
; ^= 14TeV 



-2 -1 

\ = log (X ) 



FIG. 5: Summary of the 77 — )• W~^W contribution. The lines were calculated within EPA 



approximation as described in the text with photon fluxes obtained in Ref. [4J]. Here, ^1.2 
logio(a^i,2), where x\^2 are photon longitudinal fractions with respect to parent protons. 



V. INCLUSIVE PRODUCTION OF W+W PAIRS 



For a test and for a comparison we also consider a gluon-gluon contribution to the inclusive 
cross section. We are not interested in the quark-antiquark component which is simple and 
well known. We also omit pp — )■ ttX — >■ W~^W~bbX process very important at high energy. 
In the lowest order of pQCD the inclusive cross section for the gluon-gluon fusion can be 

written as 



da^s 



dyj^dy^cPp 



wi. 



IGtt^s^ 



xig{xi, ij,p)x2g{x2, fip) \M 



gg^W+W- 



(Ai, A2, A+, A_j 



(5.1) 



The corresponding matrix elements have been discussed in the literature in detail [27| . The 
distributions in rapidity of W'^ iy+), rapidity of W~ (?/_) and transverse momentum of one 
of them pw± can be calculated in a straightforward way from Eq. (15. ip . The distribution in 
invariant mass can be then obtained by an appropriate binning. Our inclusive da/dMww 
distribution seems consistent with similar distributions presented in the past in the literature. 
The total cross section can be obtained from a simpler formula: 



a 



99 
pp- 



^w+w-x 



dxidx2g{xi, /i|) g{x2, jj^p) dgg^w+w- (■§ 



(5.2) 



Let us concentrate for a while a the elementary gg — )■ W^W cross section shown in FigjG] 



28 . We also show a vertical line at 



In this calculation we have assumed rriho = 125 GeV 

the ti threshold. The figure demonstrates a cancellation pattern between box and triangle 

contributions. We will discuss similar cancellation for the pp — )■ ppW^W~ reaction in the 
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next section. We wish to notice that CTgg^w+w- ^ cr^^^^+w- 
a potential role of photon-photon induced processes of W~^W~ 
far in the context of inclusive process. 



- ^^^ ~ 10^ pb. This shows 
production not discussed so 
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FIG. 6: The integrated elementary cross section for the gg — t- W~^W~ reaction. The solid line 
represents the coherent sum of all contributions. We show separate contributions of boxes (dashed 
line) and triangles (dotted line). 

As discussed before, in the case of exclusive scattering the J^ = contribution is the 
dominant one. In the case of inclusive process the situation is slightly different. In FigJT] 
we present the J^ = and \Jz\ = 2 components to angular distributions. The Jz = 
contribution is generally larger than the \Jz\ = 2 one. As in the exclusive case, at forward 
scattering {cos 9 = ±1) we observe the dominance of the J^ = contribution. At Vs = 500 
GeV it happens very close to cos 9 ^ ± 1. 

For completeness in Figj8]we show corresponding contributions to the rapidity distribu- 
tion of one of W^s in the pp — )■ W^W~X process. Here the Jz = contribution is larger in 
the whole range of rapidities. 

VI. RESULTS 



Before we go to the presentation of results for the pp — )■ ppW~^W reaction we wish to 
show results for the pp — )■ pp77 reaction. The latter reaction was studied experimentally in 
Ref.H- 



A. 



PP ~^ PPll 



The pp — >■ PP77 process was discussed recently in [46j. No differential distributions have 
been discussed there. The CDF Collaboration has measured photons in the interval |^(7)| < 
1.0, Erp > 2.5 GeV and with the condition of no other particles detected in -7.4 < t] < 7.4. 
They have obtained a^^ = 2.48 pb with about ^^arter of relative uncertainty. We obtain 



2.99 pb for the GJR NLO gluon distribution [47[, 2.46 pb for the MSTW08 NLO gluon 
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FIG. 7: Centre-of-mass scattering angle dependence of the hard subprocess gg — )• W~^W~ cross 
section averaged over incoming gluon polarizations. The solid line represents the coherent sum of 
all contributions The J^ = (dashed line) and the \Jz\ =2 (dotted line)) contributions are shown 
separately. 
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FIG. 8: The J^ = (dashed line) and |J^ 
pp — )• W^W^X rapidity distribution. 



2 (dotted line) contributions to the inclusive 



distribution 48| and 2.1 pb for the CT12 NLO gluon distribution [49.]. Our results very 
well agree with the CDF experimental data. In this calculation we have assumed averaged 
soft gap survival factor Sg = 0.05 and the scale of the Sudakov form factor was taken as 



/i 



77 



Cuts on the gluon transverse momenta g^ ^^^ = 0.5 GeV^ were imposed. 



In Figj9] (left panel) we show distribution of photon-photon invariant mass with exper- 
imental CDF cuts. We show results for three different gluon distributions |47H49| . We 
obtain very good description of the CDF experimental data [20], both in shape and absolute 
normalization. In the right panel we show corresponding distribution in photon transverse 
momentum. 

Finally in FigUDlwe show corresponding distribution in photon pseudorapidity in the left 
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FIG. 9: Left panel: Photon-photon invariant mass distribution. We show results for three different 
gluon distributions specified in the figure. The experimental data are taken from Ref. 20|]. Right 
panel: 



panel, again for three different gluon distributions. In the right panel we present decompo- 
sition into different pp center-of-mass photon helicity components. 
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FIG. 10: Distribution in photon pseudorapidity for three different gluon distributions (left panel) 
and the decomposition into different pp center-of-mass photon helicity components. 



Having shown that the results of the approach used in the present paper nicely describe 
the CDF experimental data 20|| we can confidentially present our predictions for the pp — )■ 
ppW^W~ reaction. 



B. pp — )- ppW^W 

Let us present now our results for the central exclusive W~^W~ pair production. In Fig. [11] 
we compare rapidity distribution of W~^ (or W~) for the electromagnetic 77 — )■ W^W~ 
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and diffractive gg — )■ W^W~ mechanisnis. The two-photon induced contribution is almost 
three orders of magnitude larger than the diffractive contribution, in which all polarization 
components for W^ and W~ have been included. For a reference, we show also inclusive cross 
section {gg — )■ W^W~ contribution only) which is roughly two more orders of magnitude 
bigger than the exclusive 77 — )• W^W~ contribution. We see, therefore, that the exclusive 
diffractive component is five orders of magnitude smaller for its inclusive counterpart. The 
diffractive contribution was calculated with the GJR NLO 47j collinear gluon distribution, 
in order to generate the off-diagonal UGDFs given by Eq. fl2.16p . This collinear PDF allows 
us to use quite small values of gluon transverse momenta (g^ ^^^ = 0.5 GeV^). 

A much smaller diffractive contribution compared to the two-photon one requires a special 
comment as it is rather exceptional. For example, it is completely opposite than for op -^ 
ppH [7|, pp — )■ ppM (e.g. light/heavy quarkonia production |9|, |lO|) or pp — )■ ppQQ [7|, |50| 
CEP processes. The standard relative suppression, present also in the latter cases, is due to 
soft gap survival probability factor {Sg ~ 0.03 for diffractive contribution versus Sg ~ 1 for 
two-photon contribution), and due to a suppression by the Sudakov form factor calculated 
at very large scales, here at fj.hard = ^vwk- The main difference compared to other cases is 
that in the diffractive case the leading contribution comes from loop diagrams while in the 
two-photon case already from tree level diagrams. 
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FIG. 11: Rapidity distribution of W bosons. The diffractive contribution is shown by the bottom 
hne while the 77 — )• W~^W~ contribution by the middle line. For comparison, we also show the 
cross section for the inclusive (55- fusion only) production case (upper line). 



In Fig. [12] we present, in addition, individual polarization components for the diffrac- 
tive mechanism, along with the unpolarized cross section. The calculation of the helicity 
contributions is performed in the pp center-of-mass frame (in which all the experimental 
studies of the exclusive production processes are usually performed). As can be seen from 
the figure, the contribution of (A+,A„) = (±1,=f1) is bigger than other contributions and 
the contribution of (A+, A_) = (±1, ±1) concentrated mostly at midrapidities. Since we use 
pp center of mass helicities there is on simple relation to the often used in a qualitative 
discussion Jz = dominance rule. Discussion of the Jz = rule would require complicated 
transformations between different reference frames and going beyond approximations made 
here. This clearly goes beyond the scope of this paper. In particular, as it is seen from 
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Fig. [12] the helicity contributions obey the following relation 

dox\'{y+) dax'\{y-) 



dy-i 



dy_ 



(6.1) 



where y± are rapidities of W^ bosons, respectively. The unpolarized cross section does not 
show up any peculiarities in ^/-dependence and is symmetric with respect to y = for both 
W^ and W~ bosons. 
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FIG. 12: Rapidity distribution of separate polarisation components to the diffractive W bosons 
production. The individual contributions are marked in the figure. 



In Fig. [13] we show distribution in W~^ {W~) transverse momentum. The distribution 
for exclusive diffractive production is much steeper than that for the electromagnetic con- 
tribution. A side remark is in order here. The diffractive contribution peaks at pt,w ~ 25 
GeV. This is somewhat smaller than for the 77 — )■ W~^W~ mechanism where the maximum 
is at pt^w ~ 40 GeV. The exclusive cross section for photon-photon contribution is at large 
transverse momenta ~ 1 TeV smaller only by one order of magnitude than the inclusive 
gg — )■ W~^W~ component. The situation could be even more favorable if New Physics 



would be at the game [21 

Fig. [13] shows distribution in the W^W' 



invariant mass which is particularly important 
for the New Physics searches at the LHC 2l|. The distribution for the diffractive component 
drops quickly with the My/w invariant mass. For reference and illustration, we show also 
distribution when the Sudakov form factors in Eq. fl2.16p is set to one. As can be seen from 
the figure, the Sudakov form factor lowers the cross section by a large factor. The damping 
is MiyvF-dependent as can be seen by comparison of the two curves. The larger Mww 
the larger the damping. We show the full result (boxes -|- triangles) and the result with 
boxes only which would be complete if the Higgs boson does not exist. At high invariant 
masses, the interference of boxes and triangles decreases the cross section. The distribution 
for the photon-photon component drops very slowly with Mww and at M\yw > 1 TeV the 
corresponding cross section is even bigger than the gg — )■ W^W~ component to inclusive 
production of W^W~ pairs. 

Finally, in Fig. [15] we show for completeness the two-dimensional distributions in rapidi- 
ties of W~^/W~ bosons in both electromagnetic and QCD mechanisms. We see a typical 
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FIG. 13: Distribution in transverse momentum of one of the W bosons. The diffractive contri- 
bution is shown by the bottom sohd hne while the 77 — )• W^W~ contribution by the middle solid 
line. The top solid line corresponds to the inclusive two-gluon initiated pp — )• W^W~X compo- 
nent. Separate contributions of boxes (dashed) and triangles (dotted) are shown in addition for 
illustrating the cancellation effect. 



pp -^ pp w"" W" 
\(s=14TeV 1 




500 



1000 1500 2000 

Mww (GeV) 



FIG. 14: Distribution in W^W~ invariant mass. We show both the QCD diffractive contribution 
and the electromagnetic 77 — )• W^W~ contribution. The result when the Sudakov form factor is 
put to one is shown for illustration of its role. The most upper curve is for the inclusive gluon- 
initiated pp — >• W^W~X component. 

correlation pattern characteristic for 2 — )• 2 subprocesses. This distribution does not show 
any specific behavior which could be used to differentiate the diffractive and the two-photon 
contributions. 
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FIG. 15: Two-dimensional distribution in rapidity of W~^ and W bosons for the diffractive 
mechanism (left panel) and two-photon mechanism (right panel). 

VII. CONCLUSIONS 

We have calculated the QCD diffractive contribution to the exclusive pp — )■ pW^W~p 
process for the first time in the literature with the full one-loop gg -^ W^W~ matrix 
element. Two mechanisms have been considered. First mechanism is a virtual (highly off- 
shell) Higgs boson production and its subsequent transformation into real W^W~ pair. 
Second mechanism relies on the formation of intermediate quark boxes, very much similar 
to ones in the exclusive two photon production mechanism. 

We have calculated corresponding amplitudes using computer program package 
FormCalc. We have made a first estimate of the cross section using amplitudes in the 
forward limit "corrected" off-forward via a simple exponential (slope dependent) extrapola- 
tion. 

In order to gain confidence to our calculations and the formalism used we consider also 
the pp — )■ PP77 process which was measured recently by the CDF Collaboration. Here the 
formalism of calculating quark box diagrams is essentially the same as for the exclusive 
production of W^W~ pairs. We have obtained very nice agreement with experimental 
diphoton invariant mass distribution. 

Having verified the formalism for diphoton production we have performed similar calcu- 
lation for W^W~ production. Differential distributions in the W^ transverse momentum, 
rapidity and W^W~ pair invariant mass have been calculated and compared with corre- 
sponding distributions for discussed in the literature 77 — )■ W^W~ mechanism. The con- 
tribution of triangles with the intermediate Higgs boson turned out to be smaller than the 
contribution of boxes taking into account recent very stringent limitations on Higgs boson 
mass from Tevatron and LHC data. We have found that, in contrast to exclusive production 
of Higgs boson or dijets, the two-photon fusion dominates over the diffractive mechanism 
for small four- momentum transfers squared in the proton lines (^1,^2) as well as in a broad 
range of iy"'"iy~-pair invariant masses, in particular, for large My/w- Estimated theoreti- 
cal uncertainties cannot disfavor this statement. The large M^y/ region is damped in the 
diffractive model via scale dependence of the Sudakov form factor. 

One could focus on the diffractive contribution by imposing lower cuts on ti and/or ^2 
using very forward detectors on both sides of the interaction point at distances of 220 m and 
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420 m as planned for future studies at ATLAS and CMS. The corresponding cross section 
is, however, expected to be extremely low. 

Compared to the previous studies in the effective field theory approach, in this work we 
have included the complete one-loop (leading order) gg — >■ W~^W~ matrix element, and have 
shown that extra box diagrams, even though they are larger than the resonant (s-channel 
Higgs) diagrams, constitute a negligibly small background for a precision study of anomalous 
couplings. 

The unique situation of the dominance of the 77 — )■ W^W~ contribution over the diffrac- 
tive one opens a possibility of independent tests of the Standard Model as far as the triple- 
boson 'jWW and quartic-boson 'j'jWW coupling is considered. It allows also for stringent 
tests of some Higgsless models as discussed already in the literature (see e.g. Ref. |2l|). 



VIII. ACKNOWLEDGMENTS 

Useful discussions with Rikard Enberg, Gunnar Ingelman, Valery Khoze, Otto Nacht- 
mann, Christophe Royon, Torbjorn Sjostrand and Marek Tasevsky are gratefully ac- 
knowledged. This study was partially supported by the MNiSW grants No. DEC- 
2011/01/N/ST2/04116 and DEC-2011/01/B/ST2/04535. Authors are grateful to the Euro- 
pean Center for Theoretical Studies in Nuclear Physics and Related Areas (ECT*, Trento, 
Italy) for warm hospitality during their stay when this work was completed. Piotr Lebiedow- 
icz is thankful to the THEP Group at Lund University (LU), Sweden, for hospitality during 
his collaboration visit at LU. 



[1] A. Schafer, O. Nachtmann and R. Schopf, Phys. Lett. B249 (1990) 331. 

[2] A. Bialas and P. V. Landshoff, Phys. Lett. B256 (1991) 540. 

[3] M. G. Albrow, T. D. Coughlin and J. R. Forshaw, Prog. Part. Nucl. Phys. 65 (2010) 149. 

[4] M. G. Albrow et al. (FP420 R&D Collaboration), JINST 4, TlOOOl (2009). 

[5] M. G. Albrow and A. Rostovtsev, arXiv:0009336 [hep-ph]. 

[6] V. A. Khoze, A. D. Martin and M. G. Ryskin, Phys. Lett. B401 (1997) 330; 

V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C14 (2000) 525; 

V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C19 (2001) 477 [Erratum-ibid. 

C20 (2001) 599]; 

V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C23 (2002) 311; 

A. B. Kaidalov, V. A. Khoze, A. D. Martin and M. G. Ryskin, Eur. Phys. J. C33 (2004) 261. 
[7] R. Maciula, R. Pasechnik and A. Szczurek, Phys. Rev. D82 (2010) 114011; 

R. Maciula, R. Pasechnik and A. Szczurek, Phys. Rev. D83 (2011) 114034. 
[8] R. Maciula, R. Pasechnik and A. Szczurek, Phys. Rev. D84 (2011) 114014. 
[9] R. S. Pasechnik, A. Szczurek and O. V. Teryaev, Phys. Rev. D78 (2008) 014007; 

R. S. Pasechnik, A. Szczurek and O. V. Teryaev, Phys. Lett. B680 (2009) 62; 

R. S. Pasechnik, A. Szczurek and O. V. Teryaev, Phys. Rev. D81 (2010) 034024; 

R. Pasechnik, A. Szczurek and O. Teryaev, Phys. Rev. D83 (2011) 074017; 

P. Lebiedowicz, R. Pasechnik and A. Szczurek, Phys. Lett. B701 (2011) 434. 
[10] L. A. Harland-Lang, V. A. Khoze, M. G. Ryskin and W. J. Stirling, Eur. Phys. J. C69 (2010) 

179. 

20 



[11] J. R. Cudell, A. Dechambre, O. F. Hernandez and I. P. Ivanov, Eur. Phys. J. C61 (2009) 369. 
[12] A. Dechambre, O. Kepka, C. Royon and R. Staszewski, Phys. Rev. D83 (2011) 054013. 
[13] V. A. Khoze, M. G. Ryskin and W. J. Stirhng, Eur. Phys. J. C48 (2006) 477. 
[14] R. Pasechnik, R. Enberg and G. Ingelman, Phys. Rev. D82 (2010) 054036; 

R. Pasechnik, R. Enberg and G. Ingelman, Phys. Lett. B695 (2011) 189. 
[15] M. G. Ryskm, A. D. Martin, V. A. Khoze, Eur. Phys. J. C60 (2009) 265. 
[16] A. Szczurek, R. S. Pasechnik and O. V. Teryaev, Phys. Rev. D75 (2007) 054021. 
[17] L. A. Harland-Lang, V. A. Khoze, M. G. Ryskin and W. J. Stirling, Eur. Phys. J. C71 (2011) 

1714. 
[18] R. Enberg and R. Pasechnik, Phys. Rev. D83 (2011) 095020. 
[19] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 99 (2007) 242002; 

T. Aaltonen et al. (CDF Cohaboration), Phys. Rev. D77 (2008) 052004; 

T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 102 (2009) 242001. 
[20] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 108 (2012) 081801. 
[21] O. Kepka and C. Royon, Phys. Rev. D78 (2008) 073005; 

E. Chapon, C. Royon and O. Kepka, Phys. Rev. D81 (2010) 074003. 
[22] N. Schul and K. Piotrzkowski, Nucl. Phys. B (Proc. Suppl.) 179-180 (2008) 289; 

T. Pierzchala and K. Piotrzkowski, Nucl. Phys. B (Proc. Suppl.) 179-180 (2008) 257. 
[23] M. Maniatis, A. v. Manteuffel and O. Nachtmann, Nucl. Phys. B (Proc. Suppl.) 179-180 

(2008) 104. 
[24] O. Nachtmann, F. Nagel, M. Pospischil and A. Utermann, Eur. Phys. J. C45 (2005) 679; Eur. 

Phys. J. C46 (2006) 93. 
[25] B. E. Cox, A. De Roeck, V. A. Khoze et al, Eur. Phys. J. C45 (2006) 401. 
[26] R. S. Gupta, Phys. Rev. D85 (2012) 014006. 
[27] E. W. N. Glover and J. J. van der Bij, Phys. Lett. B219 (1989) 488; 

Ch. Kao and D. A. Dicus, Phys. Rev. D43 (1991) 1555. 
[28] G. Aad et al. (ATLAS Collaboration), Phys. Lett. B716 (2012) 1; 

S. Chatrchyan et al. (CMS Collaboration), Phys. Lett. B716 (2012) 30. 
[29] S. Y. Choi and F. Schrempp, Phys. Lett. B272 (1991) 149; 

E. Yehudai, Phys. Rev. D44 (1991) 3434; 

G. Belanger and F. Boudjema, Phys. Lett. B288 (1992) 210. 
[30] A. Denner, S. Dittmaier and R. Schuster, Nucl. Phys. B452 (1995) 80. 
[31] K. Nakamura et al. (Particle Data Group), J. Phys. G37 (2010) 075021. 
[32] The TEVNPH Working Group for the CDF and DO Cohaborations, FERMILAB-CONF-11- 

372-E, larXiv:1108.333n [hep-ex]. 
[33] ATLAS Collaboration, ATLAS-CONF-201 1-163; CERN-PH-EP-2011-190, Phys. Rev. Lett. 

108 (2012) 111802. 
[34] CMS Collaboration, CMS PAS HIG-11-032; CMS-HIG-11-024, Phys. Lett. B710 (2012) 91. 
[35] V. I. Kuksa and R. S. Pasechnik, Int. J. Mod. Phys. A24 (2009) 5765; 

R. S. Pasechnik and V. I. Kuksa, Mod. Phys. Lett. A26 (2011) 1075. 
[36] G. Passarino, Nucl. Phys. B488 (1997) 3. 

[37] R. S. Pasechnik, O. V. Teryaev and A. Szczurek, Eur. Phys. J. C47 (2006) 429. 
[38] R. Enberg, R. Pasechnik and O. Stal, Phys. Rev. D85 (2012) 075016. 
[39] T. Hahn, Comput. Phys. Commun. 140 (2001) 418; 

T. Hahn and M. Perez- Victoria, Comput. Phys. Commun. 118 (1999) 153; 

T. Hahn, Comput. Phys. Commun. 178 (2008) 217. 

21 



[40] M. A. Kimber, A. D. Martin and M. G. Ryskin, Phys. Rev. D63 (2001) 114027. 

[41] A. D. Martin and M. G. Ryskin, Phys. Rev. D64 (2001) 094017. 

[42] T. D. Coughlin and J. R. Forshaw, JHEP 1001 (2010) 121. 

[43] P. Lebiedowicz and A. Szczurek, Phys. Rev. D81 (2010) 036003. 

[44] M. Drees and D. Zeppenfeld, Phys. Rev. D39 (1989) 2536. 

[45] M. A. Shifman, A. I. Vainshtein, M. B. Voloshin and V. I. Zakharov, Sov. J. Nucl. Phys. 30 

(1979) 711; 

E. E. Boos and G. V. Jikia, Phys. Lett. B275 (1992) 164; 

D. A. Morris, T. N. Troung and D. Zappala, Phys. Lett. B323 (1994) 421; 

H. Veltman, Z. Phys. C62 (1994) 235. 
[46] L.A. Harland-Lang, V.A. Khoze, M.G. Ryskin and W.J. Stirhng. larXiv: 1204.4803) [hep-oh]. 
[47] M. Gliick, D. Jimenez-Delgado and E. Reya, Eur. Phys. J. C53 (2008) 355; 

M. Gliick, D. Jimenez-Delgado, E. Reya and C. Schuck, Phys. Lett. B664 (2008) 133. 
[48] A.D. Martin, W.J. Stirling, R.S. Thorne and G. Watt, Eur. Phys. J. C63 (2009) 189; 

G. Watt and R.S. Thorne, arXiv: 1205.4024 [hep-ph]. 
[49] H.-L. Lai et al, Phys. Rev. D82 (2010) 074024. 
[50] R. Maciula, R. Pasechnik and A. Szczurek, Phys. Lett. B685 (2010) 165. 



22 



